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Abstract 

A model for single-domain uniaxial ferromagnetic particles with high anisotropy, 
the Ising model, is studied. Recent experimental observations have been made of 
the probability that the magnetization has not switched, -Pnot- Here an approach is 
described in which it is emphasized that a ferromagnetic particle in an unfavorable 
field is in fact a metastable system, and the switching is accomplished through 
the nucleation and subsequent growth of localized droplets. Nucleation theory is 
applied to finite systems to determine the coercivity as a function of particle size 
and to calculate -Pnot- Both of these quantities are modified by different boundary 
conditions, magnetostatic interactions, and quenched disorder. 

Key words: Magnetization - reversal. Switching Fields, Metastable Phases, 
Ferromagnets - nanoscale 



CONTACT: M.A. Novotny, SCRI, Florida State University, Tallahassee, 
FL 32306-4052; FAX 1-904-644-0098; email: novotny@scri.fsu.edu 



Preprint submitted to Elsevier Preprint 



1 February 2008 



Highly anisotropic single-domain ferromagnetic particles are candidates 
for applications in ultra-high density magnetic recording. Since such par- 
ticles do not have internal domain walls in equilibrium, wall-motion de- 
scriptions of the switching dynamics are inadequate. Recent experimental 
observations, including //-SQUID [1] and Magnetic Force Microscopy [2], 
allow measurements on individual magnetic particles of the probability 
that the magnetization has not switched at time t, Pnot(^)- Such experi- 
ments provide tests for theoretical predictions of the reversal mechanism 
for single-domain nanoscale magnetic particles. The theoretical approach 
emphasized here is that a ferromagnetic particle in an unfavorable field 
is in fact a metastable system [3] , and that the switching is accomplished 
through the nucleation and subsequent growth of localized droplets. 

To simulate the switching dynamics a kinetic Ising model with Hamilto- 
nian Ti = —JJ2{i,j) SiSj — HL'^m + DL'^m'^, is used, where Si = ±1 is the 
^-component of the local magnetization, J > for these ferromagnetic 
systems, H is the applied magnetic field times the single-spin magnetic 
moment, D is a mean- field approximation for the demagnetizing field [4], 
and a LxL square lattice {d=2) with only nearest-neighbor exchange in- 
teractions is used. The dimensionless magnetization is m = L~'^J2iSij 
where the sum runs over all L'^ spins. The dynamic is to choose a spin 
at random, and then flip it with the Metropolis or with the Glauber 
spin-flip probability [3] . Time is measured in Monte Carlo steps per spin, 
MCSS. The simulation starts with all spins +1 and a negative field H, 
then Pnot(^) is measured. Fixing a waiting time, t„, the switching field, 
Hsw, defined as the field at which Pnot(^w)=|, is also measured. 

The simulations are performed at temperatures well below the critical 
temperature for the pure system, Tc~2.269J. There are at least four rel- 
evant length scales in the problem [3]: the lattice spacing (which is set 
equal to unity); the system size L; the radius of a critical droplet. Re, and 
the typical distance between droplets that are critical or supercritical, 
Ro- This paper concentrates on three regimes [3,5]: the coexistence (CE) 
regime where L < R^; the single-droplet (SD) regime where Rc < L < Rg] 
and the multi-droplet (MD) regime where Ro < L. Figure 1 presents ifgw 
as a function of L for various situations. The curves with bond dilution 
were simulated with the Glauber dynamic, while the Metropolis dynamic 
was used in the other curves. On this scale results interchanging the dy- 
namics should be indistinguishable. The maxima in ifgw are near the 
cross-over between the GE and SD regimes. The shapes of these curves 
for periodic [5] and other [6] boundary conditions have been predicted by 
detailed droplet-theory calculations. Results for quenched disorder have 
been presented recently [7]. 

In the SD regime [5] Pnot(^) = exp(— i/r), where r is the average lifetime 
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of the metastable state. However, in the MD regime [5] 



PNOT{t) w - 1 - erf 



t-r 




(1) 



where for D=0 and periodic boundary conditions At has been obtained 
exphcitly [5]. Figure 2 shows Pnot(^) in the MD regime. The shape of 
-Pnot(^) is not changed for finite D or for randomness in the field. 

In summary, magnetization switching in a finite system is comphcated, 
even for the simplest highly anisotropic ferromagnetic model. There are 
always at least four length scales in the problem, and the switching mech- 
anism is different in regimes with different relationships between these 
length scales. Even with no demagnetizing fields highly anisotropic single- 
domain magnetic particles exhibit a maximum in the switching field as a 
function of system size. The form for Pnot is different in different regimes: 
in the SD regime it is an exponential while in the MD regime it is well 
approximated by an error function. 
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Figure Captions 

Fig. 1. The switching field with ^=0, tw=30000 MCSS, and temperature 
T=1.3J. SoUd fines witfi x are for periodic boundary conditions and 
random bond dilution. Otfier data are for a circular system with free 
boundary conditions (o) , and for a square system with periodic boundary 
conditions in one direction and free boundary conditions in the other 
direction (□). 

Fig. 2. The probability of not switching, Pnot(^)) in the multi-droplet 
(MD) regime for periodic boundary conditions; L=100, T=0.8Tc~1.815J, 
using Glauber dynamics. The curves are from 1000 escapes from the 
metastable state. Curve (a) for //=-0.34725 J and ^=0; (b) for //=-0.34725 J 
and D—0.05J, and (c) for D—0 and a uniformly distributed random field 
centered at H=-0M725J with a width of 0.34725J. This figure should 
be compared with Fig. 3c of [5] . which is for the Metropolis dynamic, but 
has the same values for L, T, and H. 
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